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A  DEFINITION  OF  STABLE  INELASTIC  MATERIAL* 
D.  C.  Drucker**  (Brown  University) 


Abstract 

The  definitions  given  previously  for  work-hardening  and  perfect 
plasticity  are  broadened  to  cover  viscous  effects.  As  before,  the  system 
con/sidered  includes  both  the  body  and  the  time-dependent  set  of  forces 
acting  upon  it#  An  external  agency  is  supposed  to  apply  an  additional 
set  of  forces  to  the  already  loaded  body#  It  is  now  postulated  that 
the  work  done  by  the  external  agency  on  the  additional  displacements  it 
produces  is  positive*  Some  of  the  restrictions  thus  imposed  on  per¬ 
missible  stress-strain  relations  are  explored*  Especial  attention  is 
paid  to  simple  laws  of  creep*  Uniqueness  of  solution  also  is  studied* 


*  The  results  presented  in  this  paper  were  obtained  in  the  course  of 
research  sponsored  by  the  Office  of  ll&val  Research  under  Contract 
lkmr«562{20}  with  Brown  University* 

Cbadman*  Division  of  Ihgteeertng* 
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Introduction 

In  the  description  of  material,  terms  such  as  elastic,  viscous, 
and  plastic  are  eiqployed  separately  or  in  combination  as  in  elastic-plastic, 
visco-elastic,  and  visco-plastic.  Their  meaning  is  quite  clear  in  a 
qualitative  sense.  The  difficulty  lies  in  obtaining  a  broad  and  physically 
valid  definition  of  terms  which  will  lead  to  a  useful  mathematical 
formulation  of  general  stress-strain  relations.  Any  restrictions  imposed 
limit  the  class  of  materials  which  are  described  and  so  are  arbitrary 
to  a  considerable  extents  In  this  paper  and  its  predecessors  (1)  (2)* 
the  restriction  is  to  stability  of  material  in  the  strict  sense. 

An  elastic  body  is  reversible  and  non-dis  sipat ive  and  time 
independent  under  isothermal  conditions.  All  work  done  on  the  body  is 
stored  as  strain  energy  and  can  be  recovered  on  unloading.  The  stress- 
strain  curve  in  simple  tension  may  have  a  wide  variety  of  shapes  including 
the  unstable  portions  shown  in  Figure  1,  If  stable  elastic  material  is 
specified  the  drop  in  stress  with  increase  in  strain  must  be  ruled  out 
explicitly. 

Plasticity  denotes  irreversibility  and  permanent  or  residual 
strain  upon  unloading.  In  the  narrow  but  convenient  terminology  of  the 
mathematical  theory  of  plasticity,  time  independence  of  material  properties 
is  understood  despite  the  frequent  us#  of  terras  such  a/3  strain  rate  and 
velocity.  Time  always  enters  homogeneously  in  all  equations  of  statics, 

A  given  loading  path  produces  the  same  deformations  regardless  of  the 
actual  time  history  of  the  loading.  Figure  2  gives  possible  simple  tension 
curves  for  a  material  combining  linear  elasticity  and  work-hardening 
plasticity.  The  phenomenon  of  the  upper  yield  point  in  mild  steel  1* 

^waw>«wr.>wiBgw»aiiwraiw|l|(i  nwrurtyn  m  jmwtmiii  r  iW'Pw.Mii'iBWL  jojulw  i  i(riji.liiiMnaM,ii r  t nji wwn 

Numbers  in  parentheses  refer  to  tee  Bibliography  at  the  end  of  the 
paper# 
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a  convincing  demonstration  of  the  fact  that  instability  of  material  in  the 
mechanical  sense  is  encountered  in  the  physical  world*  A  gradually  falling 
stress-strain  curve  likewise  is  a  physical  possibility  although  not  found 
in  structural  metals  when  actual  rather  than  nominal  stress  is  plotted* 

If  stable  elastic-plastic  material  is  specified,  therefore,  the  drop  in 
stress  must  be  ruled  out  explicitly* 

Time  effects  are  present  to  some  extent  in  all  real  materials* 
When  the  stress-strain  curve  is  affected  significantly  by  the  rate  of 
loading  or  more  generally  by  the  time  history  of  the  loading,  where  time 
enters  in  an  essential  manner,  the  viscosity  or  the  modifier  visco- 
is  used*  The  idealization  known  as  linear  visco-elasticity  may  be 
visualized  in  terms  ©f  combinations  of  linearly  elastic  springs  and 
Newtonian  dashpots,  Figure  3*  Non-linearity  of  the  springs  and  the 
dashpots  poses  no  conceptual  difficulty  nor  does  the  inclusion  of  elastic- 
plastic,  "spring51  elements*  Falling  stress-strain  curves  are  easily  pro¬ 
duced  by  a  simple  time  history  in  the  most  eleoientaiy  Kelvin  or  Maxwell 
models*  Instability  now  has  a  more  obscure  meaning  and  certainly  must 
involve  time*  For  example,  a  falling  curve  on  a  plot  of  stress  versus 
strain-rate  for  a  Maxwell  type  material  would  denote  viscous  instability. 
Figure  3* 

Previous  Postulate  for  Elastic  and  Plastic  Materials* 

A  fundamental  definition  has  been  advanced  previously  for  a 
stable  elastic  or  work-hardening  tii^independent  material  (1)  (2)* 
Consider  a  body  of  such  material  at  rest  acted  u pm  by  a  set  of  boundary 
tractions  and  body  forces  Figure  An  external  agency  is 
to  aid  slowly  a  set  of  surface  tractions  A %  and  body  forces 


Honr-562(20)/2 


and  equilibrium  is  maintained*  It  was  postulated  that  the  work  done  by 
the  external  agency  during  the  application  of  force  must  be  positive  and 
over  a  cycle  of  application  and  removal,  positive  or  zero,  zero  only  when 
purely  elastic  changes  take  place*  In  such  a  cycle,  work  cannot  be 
extracted  from  the  system  of  the  body  and  the  set  of  forces  acting 
upon  it*  Materials  following  this  postulated  behavior  are  stable*  A 
falling  stress-strain  curve  in  simple  tension  is  ruled  out  as  is  far  more 
complicated  behavior  which  normally  might  not  be  considered  as  unstable* 
Despite  the  severity  of  the  definition,  the  consequences  do  seem  in  accord 
with  all  existing  intuitive  concepts* 

The  implications  or  consequences  of  the  fundamental  definition 
are  remarkably  strong  for  the  mathematical  theory  of  plasticity*  The 
yield  or  subsequent  loading  surface  f(o^j)  =  which  separates  states 
of  stress  which  can  be  Joined  by  a  purely  elastic  path  from  those  requiring 
elastic-plastic  action.  Figure  5,  must  be  convex*  Furthermore,  as  shown, 
the  plastic  strain  increment  or  strain  "rate”  must  be  normal  to  the  surface 
at  a  smooth  point  =  Xdf/dcry  or  lie  between  the  outward  normals  to 
adjacent  points  at  a  comer* 

Both  convexity  and  normality  had  been  postulated  long  before 
(3)  (4)  and  it  is  certainly  valid  to  question  the  value  of  replacing 
one  set  of  postulates  by  another*  One  answer  is  that  it  is  desirable  to 
have  the  smallest  number  of  postulates*  A  better  statement,  perhaps,  is 
that  fundamental  definitions  have  many  specific  implications  beyond  those 
of  immediate  use#  Direct  experimental  study  and  theoretical  consideration 
of  a  basic  postulate  form  the  basis  of  unlimited  future 
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Coulomb  type  do  not  follow  the  postulate  (5).  A  sub-class  of  materials, 
therefore,  is  specified  and  the  problem  reduces  to  how  well  particular 
solids  are  described* 


An  Extended  Postulate  or  Definition 

The  absence  of  time  effects  in  the  mathematical  theory  of 
elasticity  and  plasticity  is  not  matched  in  the  physical  world  (6)* 

Bodies  under  constant  load  will  creep,  bodies  under  constant  deformation 
will  relax*  For  many  structural  metals  at  working  temperatures  the 
time  effects  are  small  enough  to  be  ignored  in  practical  applications 
and  the  previous  postulate  is  adequate*  However,  all  metals  at  elevated 
temperatures  and  most  materials  at  room  temperature  display  appreciable 
time  effects.  A  modified  and  broadened  postulate,  therefore,  is  required* 
When  time  or  viscous  effects  are  appreciable,  it  is  no  longer 
true  that  the  work  done  by  an  external  agency  need  be  positive*  Spring 
elements  as  in  Figure  3  generally  will  contain  energy  which  can  be  extracted 
from  the  system.  A  body  under  constant  load  usually  does  not  maintain 
an  unchanged  configuration.  Application  of  additional  force  to  an  already 
loaded  body  by  an  external  agency  is  accompanied  by  a  change  in  the  dis¬ 
placement  history,  partly  due  to  the  existing  force  system  and  partly 
produced  by  the  external  agency*  Figure  4  may  now  be  reinterpreted. 

Figure  4a  shows  the  surface  tractions  %  and  body  forces  F^  which 
are  functions  of  time*  The  displacements  u^,  the  stresses  cr^,  and 
the  strains  are  likewise  functions  of  time*  Figure  4b  shows  the 

existing  system  plus  the  external  agency  and  the  combined,  effect*  The 
added  surface  tractions  AT^  and  the  added  body  forces  tW^  also  vary  I 
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superposed  on  the  existing  set  of  forces*  Except  for  the  very  special 
fully  linear  case,  it  will  not  be  true  that  AT^  and  AF^  acting  alone 
(T^  =  0,  =  0)  would  cause  Au^,  Ao^j,  ^  *^s  ^i  +  an(* 

Fi  +  AF^  which  produce  +  Au^r  +  ^aij»  €il  +^®i^»  a^*  1W^LC^ 

are  functions  of  time# 

A  fundamental  definition  of  a  stable  inelastic  ( elastic— vxsco« 

» 

plastic)  material  can  now  be  stated; 

The  work  done  by  the  external  agency  on  the  change 
in  displacements  it  produces  must  be  positive  or  zero» 

It  is  seen  immediately  that  the  previous  postulate  is  included  in  this  , 
broader  one  because  no  change  in  displacement  takes  place  in  an  elastic¬ 
workhardening  material  at  constant  load. 

The  definition  may  be  put  in  mathematical  form  as 


Jt=0 


w 

|  J  A TjAUjdA  +  J  a^AujdV 


dt  >  0 


[1] 


where  the  dot  indicates  derivative  with  respect  to  time,  and  t  =  0  is 


the  instant  of  time  at  which  the  external  agency  begins  to  apply  force# 

The  summation  convention  is  followed  for  repeated  subscripts  in  the  same 

•  *  *  * 

term,  for  example  AT^  Au^  =  ATX  Au^,  +  ATy  AUy  +  ATZ  AUZ  and,  for  later 

y  6y  •  ~z  ~s  xy  *  xy  yz  *  yz  zx  -  Zx 


<>ij  Ho  53  CTx  ex  +  av  *y  +  a*  «*  +  rw  +  + 


'X  wx 

in  the  usual  notation. 


It  is  sc-iuetimes  more  instructive  to  think  of  two  alternative 
paths  of  loading  T^,  Fp^t  and  T^,  which  diverge  at  t  s  0, 

Figure  6*  Each  would  produce  its  own  u^f  s  The  definition  thsn 
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Jt^O 


(l)ir*(2)  •(!> 


1  [u 


-  Uj 


']  dA 


■f 


Jv[Fi2)  "  ["i2) 


dt  >  0 


which  is  really  identical  to  1 1] «  It  has  the  advantage  of  focussing 
attention  on  the  meaning  of  the  incremental  quantities  and  the  type  of 
stability  which  is  postulated* 


Stability  in  the  Large  and  Small  -  Permissibility  of  Path 

Stability  in  the  large  and  stability  in  the  small  correspond 
respectively  to  tc  unrestricted  and  to  tc  limited  to  the  immediate 
neighborhood  of  t  =  0*  As  has  been  stated  the  requirement  of  stability 
is  an  arbitrary  one  in  a  sense*  The  requirement  that  the  tensile  stress- 
strain  curve  for  an  elastic-plastic  body  be  a  continuously  rising  one 
(  complete  stability)  is  not  of  real  consequence  when  the  stress  point 
is  on  a  portion  of  the  curve  which  rises  and  small  (infinitesimal)  changes 
in  stress  are  examined  (stability  in  the  small)* 

There  is  a  further  choice  involved  within  each  classification 
of  stability*  Any  path  labelled  (1)  in  Figure  6  may  be  permitted  or 
instead  T^  and  Fp^  may  be  restricted  to  constant  values,  their 
values  at  t  *  0.  It  must  be  kept  in  mind  always  that  inelastic 
materials  are  irreversible  and  not  all  paths  are  permissible  paths  in 
general#  For  example,  at  t  ~  t0  the  system  will  be  at  A  if  it  follows 
path  (X)  and  at  B  if  it  follows  path  (2)®  With  few  exceptions;  a  straight 
Una  path  from  A  to  8,  t  a(T^  -  f|^)#  »  f|^  + 

£«i25  -  Fx’ht  with  a  going  from  0  to  1,  will  not  change  the 
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the  straightline  will  exist  which  is  capable  of  giving  the  sane  final  state 
B*  In  fact  for  the  most  general  materials  there  is  no  available  path  what¬ 
soever  from  (l)  to  (2);  path  (2)  is  the  only  way  of  producing  the  state 
of  the  material  at  B* 

In  view  of  the  irreversibility  there  is  no  clear  decision  to 
be  made  on  how  much  stability  to  require*  Stability  in  the  small  with 
path  (1)  restricted  to  the  point  corresponding  to  t  =  0  would  seem 
essential*  The  generalized  definition  without  restriction,  or  stability 
of  material  in  the  complete  sense,  seems  desirable.  Some  of  its  consequences 
will  be  explored  in  what  follows. 


Stress-strain  Relations 

The  consequences  of  the  basic  postulate  with  respect  to  stress- 
strain  relations  are  strongest  and  seen  most  easily  for  a  homogeneous  system 
whose  mass  or  kinetic  energy  can  be  neglected*  Surface  tractions,  body 

forces,  and  velocities  in  Inequality  [2]  then  can  be  replaced  by  stresses 

« 

and  strain  rates  through  the  theorem  of  virtual  work.  For  any 

continuous  u* 

<5  «ij  "  ™ 


f  .  . 

]ATiuidA  + 


'J 


*  f 

h  \  & 


where  T^,  FJ,  are  any  set  of  equilibrium  values  and  V  Cy  are 

any  compatible  set  (7)*  There  need  be  no  relation  between  the  equilibrium 


and  the  compatible  sets*  Substitution  of 


(2)  _  .(1) 


for  T 


*  *  (2)  *(1) 

vr  J  -  u. 


for  etc.,  in  the  virtual  work  expression  [3]  thus  is  permissible 

(?)  ( 1 )  (2)  (1.) 

although  in  general  Tj' *  -  T}  ,  -  F*  '  alone  will  not  produce 

*(2)  '(I) 


Integration  over  the  volume  can  be  omitted  for  a  homogeneous 
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state  of  stress  and  strain  with  the  result 


ft, 


Jt=G 


[or 


(2) 

U 


(2) 

U 


>  0 


[*] 


In  the  neighborhood  of  t  =  0  the  stresses  and  strain  rates 
can  be  expanded  in  a  power  series  for  each  of  the  paths  (1)  and  (2),  the 
subscript  zero  denoting  the  value  at  t  =  0. 


*lj  =  aul +;iJ]ot+ 


2/?  +  ... 


e 


o 

•  •• 


*  •*  2, 

ij  =  6ijt+  *i310t+  4  +  •” 


t  '2 
2 


[5] 


The  two  paths  are  identical  up  to  t  «  0  so  that  stress  and  strain  are 

the  same  at  t  «  0,  ofV]  an<*  tfVi  ®  e*4^]  •  derlva- 

id  Jo  ij  Jo  ij  Jo  o 

tives  of  these  quantities  need  not  be  the  same  so  that  strain  rates  and 
stress  rates  may  differ.  For  tc  very  close  to  t  =  0,  Inequality  [4] 
gives 

*°ij  VJ id  ij  1.  V2 


r  *(2)  •(!),  fv<2)  •;(!),  +  i  f*'(2)  V(1)l  r*(2)  *C1>,  1  1 3/  -> 

laij  3o  [eij  ■  eij  ]o  2  1  ij  ij  JoUi3  "eij  Jo  '  c/  3 


+  • 


>  0 


[6] 


As  tc  may  be  made  vanishingly  small,  it  is  necessary  that 
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a?X0,» 


If,  however,  e 


!(2)  !(1) 


ij  =  eiJ 


then  it  is  necessary  that 


r„‘<2)  -a(1)i  rV(2>  -  V<1)l  >  0 

lo.,  i0^eij  6  in  -  0 


ij 


[81 


This  process  may  be  continued  to  higher  and  higher  time  derivatives  of 

strain  and  of  stress  as  well  if  etc.  The  general  require* 

*3  ~  ij 

went  in  the  sense  of  [7]  and  [8]  is 


[ 


(2) 
8  CT1.1 

3tm 


m  (1) 

s  alj 

3tm 


.n  (2) 

Ull 

dtn 


n  (1) 
d  t 


-iial]  >  o 

3tn  Jo  - 


19] 


with  m  and  n  each  >  1*  Inequality  [  9]  insures  stability  in  the 
small  in  the  broad  sense. 

(i]  f\\ 

If  a.  is  restricted  to  remain  constant  at  o' /i  ,  as  time 
ij  ij  Jo 

goes  on, for  the  more  restricted  definition  of  stability  in  the  small,  [7] 
and  [  8  ]  are  replaced  by 


(2)  »(2)  «(l) 


°i j  Jo  W 3  ‘  ®ij  Jo  -  u 


[10] 


•  (2)  •  6  (2)  •  •  (1), 


[11] 


For  time  independent  materials  the  (1)  terms  disappear  completely  leaving 

•  • 

the  earlier  result  ^  >  0. 

Return  now  to  stability  in  the  large  and  to  time  dependent 


materials*  If  the  path  QA  of  Figure  6  is  traversed  very  rapidly  and 
then  the  system  remains  at  A,  and  a  similar  time  history  holds  for 
path  OB,  the  fundamental  postulate  in  the  form  [4]  requires 
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*  B  *  A 

In  general,  the  strain  rates  e^j  and  will  depend,  of  course, 

upon  this  specially  selected  time  history#  Here  the  distinction  between 
the  broad  and  narrow  derinition  of  stability  does  not  arise  because  the 
path  OAB  traversed  with  extreme  rapidity  gives  the  same  state  point  B 
as  path  OB,  for  materials  which  need  time  to  deform# 

Many  special  paths  may  be  considered  but  their  treatment  will 
be  deferred#  Specific  results  for  plasticity,  linear  visco-elasticity, 
and  creep  will  be  discussed  next# 


Consequences  for  Plasticity 

As  has  been  mentioned,  the  previous  postulate  for  time- 
independent  materials,  with  held  constant,  has  strong  implications 

in  plasticity.  All  yield  and  subsequent  loading  surfaces 
history  of  loading)  *  k  must  be  convex#  The  plastic  strain;  rate  vector 

*p 

e .  .  must  be  normal  to  the  surface  in  stress  space.  Figure  5*  in  the 

9  Af*  * 

extended  sense*  At  a  smooth  point  of  f  =  k~,  for  -xs—  <j  >  0  or 

damv,  mn 

I  ***** 

an  outward  pointing  a  which  produces  plastic  deformation# 


=  q  _§£_  0 

iJ  SOy  3°mn 


[13] 


where  G  may  be  a  function  of  stress,  strain,  and  the  prior  history 
of  loading#  It  may  also  depend  upon  the  stress  rate  o^#  Time  is 
not  truly  relevant  so  that  G  must  be  homogeneous  of  order  zero  in  a, 

,»  i»  9 


,  G  =  f"  [1  + 


(a  amJ3 
Wmnuran/ 


]  • 


The  broader  definition  in  this  paper,  of  stability  in  the  small, 
restricts  G  still  more#  It  is  permissible  to  deal  with  plastic  strain 
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rates  only  and  substitution  in  Inequality  [7]  gives 


t 


*(2)  af 

ff  W  «MKH» 

ij 


•a;.  . 

*13 ]  [ 


*(2) 

a 

inn 


,(1)  3f  *(1) 


do 


isn 


]  >  0 


an 


[14] 


The  magnitude  of  the  stress  rate  is  arbitrary  and  can  be  chosen  to  give 
a  constant  value  of  the  component  normal  to  the  yield  or 

loading  surface#  Suppose  this  is  done  for  each  pair  (1)  and  (2) 
so  that  the  left  bracket  of  [14]  is  zero#  The  bracketed  expression 


containing  G^  and  G^ 


then  must  likewise  be  zero#  If  not,  a  very 

•(2) 

small  increase  or  decrease  in  the  magnitude  of  or  the  opposite 

for  crj^  will  violate  Inequality  [  14]  •  Therefore,  =  G^  or 

G  is  independent  of  o^#  The  stress-strain  relation  [  133  thus  is 
restricted  to  linearity  in  the  increments  or  rates  of  stress  and  strain# 
All  incremental  stress-strain  relations  now  in  use  are  linear  or  are 


combinations  of  linear  forms  as  a  matter  of  convenience#  Although  no 


pro-conceived  notions  are  upset,  it  is  of  great  interest  and  possibly 
of  value  to  know  that  a  single  postulate  or  definition  leads  inevitably 
to  all  the  main  features  of  plastic  stress-strain  relations# 

Stability  in  the  small  alone  was  employed  in  this  application 
of  the  postulate,  but  stability  in  the  large  is  needed  for  convexity. 

As  stated  before,  there  is  always  the  possibility  that  stability  in  the 
large  with  the  new  postulate  may  be  more  restrictive  than  is  reasonable 
on  physical  grounds  for  all  paths  of  loading  in  all  materials* 


Consequences  for  Linear  Visco-elasticity 

A  linear  visco-elastic  material  may  be  defined  in  terms  of 
nodels.  Figure  %  with  linear  springs  and  dashpots  in  combination  (8)# 

An  equivalent  definition  may  be  given  by  e m  or  wore  differential  equations 
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In  vhleh  tiirie  derivatives  of  any  order  may  appear  but  stress  and  strain 

enter  linearly*  The  difference  between  two  solutions  of  such  equations 

(1) 

i m  a  solution  as  well*  Therefore,  if  stress  state  and  history  a  ' 

•(1)  (2)  •(?)  * 

produce  e..  and  stress  state  and  history  o  produce  , 

(2)  (i\  *(2)  *(1) 

cr:*  -  a)1,'  will  give  c .  ,  -  s.  .#  A  stable  linear  viscoelastic  material 

ij  ij  ij  ij 

according  to  definition  as  expressed  by  the  basic  Inequality  [4]  must 
obey 


t/t=0 


aij  6ij dt  - 0 


[15] 


This  is  no  surprise  at  all  as  it  is  simply  a  statement  that 
work  must  be  done  on  an  unstressed  and  unstrained  viscoelastic  material 
in  order  to  deform  it*  Ary  combination  of  linear  springs  and  dashpots 
thus  provides  a  stable  material#  Although  the  result  is  not  unexpected, 
it  is  comforting  to  know  that  the  implications  of  the  broadened  definition 
are  in  accord  with  physical  intuition# 


Consequences  for  Creep 

. . — rrm  sa  i1  > 

Non-linearity  is  the  rule  rather  than  the  exception  in  inelastic 
raate, rials#  The  creep  of  metals  in  particular  is  markedly  non-linear# 
Consider  now  a  Maxwell  type  material,  Figure  35which  deforms  indefinitely 
under  constant  load*  The  elastic  response  can  be  ignored  in  discussing 
the  viscous*  Suppose  that,  for  convenience  and  for  reasonable  agree¬ 
ment  with  actual  behavior,  the  viscous  strain  rate  $Y.  is  a  function 

,  *v 

of  stress  alone*  As  the  existing  state  of  stress  determines  the  strain 
rate,  the  stress  rate  has  no  immediate  influence!  lt*ij  “eiJ  J 

in  Inequality  [7]  is  automatically  aero  at  any  instant  of  time  when 
two  paths  of  loading  diverge#  Inequality  [8]  for  stability  in  the  snail, 
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stability  in  the  large*  gives  far  more  spectacular  and  useful  results* 

Referring  to  Inequality  [12],  the  restriction  on  the  path  or  time 

B  A  *3  *A 

history  of  loading  is  lifted*  (c^  •  cr^K  s  ^  -  £^j)  >  0  applies  to 

» 

any  two  stress  states  at  any  tine  because  c  is  supposed  to  depend  upon 

A 

a  alone  and  so  is  not  path  dependent#  In  particular,  may  be  chosen 

as  the  stress  free  state  and  the  result  is 


JO  *Q  •  ■Q 

°ij 


[16] 


or  the  rate  of  doing  work  is  always  positive*  As  the  strain  rate  is 
viscous,  and,  therefore,  dissipative,  this  too  is  an  obvious  result  where¬ 
as  [12]  is  not* 

The  restrictions  on  stress-strain  relations  in  creep  imposed 

by  [12]  and  [l6]  are  strong  and  may  be  explored  for  each  special 

case.  In  pictorial  terms  as  shown  in  Figure®  7,  8,  the  strain  rate  vector 

3  ®  p 

makes  an  acute  angle  with  the  stress  vector,  a4  ^  e^4  >  0;  the  increment 

— *i  i»j  *” 

in  strain  rate  for  any  stress  path  is  in  the  direction  of  the  change  in 

#3  A  *  x 

stress  (ay  >  °* 

Suppose,  for  example;  that  it  is  possible  to  write 


*  _  dp 


[17] 


where  q>  is  a  function  of  the  existing  state  of  stress  only*  The  surface 
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sense  |  as  at  a  comer  of  such  a  surface  likewise  inquires  convexity  if 
[12]  is  to  be  satisfied,  Figure  8b® 

k  common  assumption  for  creep  in  simple  tension  (not  necessarily 
in  accord  with  the  actual  physical  behavior)  is  that  the  strain  rate  is 
proportional  to  some  power  n  of  th©  stress  (9)®  A  generalization  of 
this  hypothesis  is  to  choose  <p  in  [l?J  as  a  homogeneous  function  of 
stress  of  degree  n  +  1*  Substituting  in  [l6]t  the  rate  of  doing  work 
or  rate  of  dissipation  is  from  Euler *s  theorem  on  homogeneous  functions 

w(oij)  5  °ij  =  aij  35^  =  (n+1)<P  [18] 


Under  these  conditions  (p  =  constant  is  a  surface  in  stress  space  of 

constant  rate  of  dissipation  of  energy.  Therefore  W  ~  constant  defines 

a  convex  surface  in  stress  space® 

The  additional  assumption  of  isotropy  requires  W  to  be  a 

function  of  the  stress  invariants,  for  example,  the  Siam  of  the  principal 

1  1 

stresses  and  the  invariants  J2  58  2  sij  sij*  "  3  sij  sjk  ski 

the  stress  deviation  tensor  s^,  (10).  Considerable  discussion  of  convex 
homogeneous  functions  of  stress  with  these  invariants  has  been  given  in 
papers  on  plasticity  (10) (11)  and  need  not  be  repeated  here®  As  an 
example,  the  requirement  of  convexity  restricts  the  constant  c  in 
(n+l)(p=  W  =  b(J^  -  cJ3)  to  no  more  than  2.25* 


Uniqueness  of  Solution  to  Equilibrium  BrouMary  Val^e  Problems  (12) (13) 
Consider  a  body  at  time  t  «  0  whose  state  at  each  point  can 
be  specified  completely  as,  for  example.,  the  initial  unstressed  and 
unstrained  condition®  Suppose  surface  tractions  are  specified  on 
the  portion  of  the  surface  jLj»,  displacements  %  on  the  remaining 
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portion  and  body  forces  throughout  the  volume  V,  each  given 

as  a  function  of  time*  The  question  of  uniqueness  concerns  the  possibility 
of  two  (or  more)  solutions  (l)  and  (?)  for  the  interior  stresses*  strains* 
and  displacements  or  any  of  their  time  derivatives. 

Assume  that  two  equilibrium  solutions  do  exist  which  satisfy  the 
boundary  conditions*  Substitute  their  difference  in  the  equation  of  virtual 
work  [3]  just  as  in  the  description  following  Equation  [3]*  The  general 
result  is  that  the  left  hand  side  is  always  zero  and  for  any  order  of  time 
derivative  of  stress  and  strain 


8 

atm 


(2) 

6ij 


[19] 


which  is  reminiscent  of  Inequality  ’[9]  except  that  here  the  notation 
is  meant  to  include  the  zero  order  as  well*  i*e* 

,  (2)  (l),r  (2)  (IX  r  (2)  (1)  *(2)  *0-),  0+„ 

t°ij  Kj  '  eij  3  ’  etc* 

„  .  ...  ,  am  (2)  3°  _U)  „ 

Uniqueness  in  the  sense  of  — -  a.  .  =  —  e. or 

„  ,  ,  3tra  lJ  3ta 

An  *2)  (x) 

-2—  ej  J  ~  -  e  4  4  is  assured  if  the  integrand  of  [193  is  positive. 

0tn  10  d  tn  13 

Just  as  in  the  special  case  of  plasticity*  a  discussion  of  uniqueness  in 
such  great  generality  cannot  be  more  specific  (13)*  l^ck  of  uniqueness 
which  may  be  present  in  the  stress-strain  relation  appears  here  as  well® 

With  this  understanding  the  term  uniqueness  will  be  used  without  qualification 
in  the  following* 

Reference  to  Inequality  [7]  shows  that  stability  in  the  small 

*  *  • 
guarantees  uniqueness  of  a  or  e  at  each  instant  of  time*  If  c.  . 

,  ij  00J 

is  dependent  on  a,  ,  and  net  o4  4t 


p 
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of  importance  in  the  study  of  dynamic  problems  in  which  mass  or  inertia 
cannot  be  neglected.  If  needed,  this  process  can  be  continued  to  any  term 
of  Inequality  [9]*  The  initial  state  and  the  complete  history  of  loading, 
therefore,  determine  the  stress  and  the  strain  uniquely  throughout  any 
body  composed  of  inelastic  material  obeying  the  definition  of  stability 
in  the  small. 

It  might  be  ejected  that  stability  in  the  large  guarantees 
uniqueness  of  stress  independent  of  the  path  of  loading*  This  is  true 
for  any  material  undergoing  the  rapid  loading  type  of  path  described  for 
Inequality  [12]®  It  also  is  true  for  those  stable  materials  following 
any  creep  law  in  which  the  strain  rate  depends  upon  the  stress  only  because 
they  too  most  obey  Inequality  [  12  )•  It  is  not  true  for  most  materials* 

Certainly  the  fundamental  definition  does  require  that  if  a 
path  exists  from  A  to  B,  Figure  6,  the  work  done  by  the  external  agency 


on  the  change  in  displacements  it  produces  must  be  positive  or  zero, 

. *B  *  «  a 

J  ^crij-cij)(eij- eij)  ^  However,  in  general  the  integral  bears  little 

A  B  A  *  B  " 

or  no  relevance  to  the  virtual  work  term  whose  positive¬ 

ness  assures  uniqueness.  Furthermore,  no  path  AB  is  a  permissible  one 


in  the  general  case. 


Hyrianvi  qq 

All  of  the  emphasis  so  far  has  been  on  equilibrium  problems*  Hass 
and  kinetic  energy  have  been  ignored,  although  velocity  and  acceleration 
terms  are  finite*  As  will  be  seen,  the  inclusion  of  dynamic  terms  does  not 
affect  the  conclusions  to  be  drawn  and,  in  fact,  makes  less  demands  on  the 
stress-strain  relations*  Homogeneous  states  of  stress  and  strain  are  not 
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more  inclusive 


where  the  second  term  is  the  value  of  the  volume  integral  at  tc  minus 

the  value  at  the  time  of  divergence  of  paths,  t  =  0.  If  the  external 

•• 

agency  applies  load  slowly,  the  acceleration  u^  is  the  same  for  the 
two  paths  (1)  and  (2),  Figure  6*  Necessarily,  in  all  cases  the  velocity 
ii£  }  throughout  the  volume  is  the  same  as  u  at  t  =  0  and  the 

second  term  of  [20]  is  positive  or  zero  at  all  times  tc*  The  Inequality 
[4]  is  much  more  significant ,  therefore,  than  [20]* 

Conclusions 

A  single  postulate  or  definition  is  advanced  for  all  stable  in¬ 
elastic  materials*  Its  implicatic&is  are  far  reaching*  In  plasticity 
theory,  which  is  a  particularly  well-explored  subject,  the  concepts  of 
convexity,  normality,  and  incremental  linearity  all  are  consequences  of 

the  definition.  An  interesting  res^ilt  also  is  obtained  for  a  generalization 

*  n 

of  the  simple  tension  creep  law  e  =  Kd  •  The  surface  of  constant  rate  of 
dissipation  for  the  special  case  is  convex  and  e  is  nozmal  to  it  just  as 
for  the  yield  surface  in  plasticity*  However,  non-linear  viscous  behavior 
coupled  with  elastic  and  plastic  action  is  an  almost  unknown  field*  liich 
more  must  be  done  before  the  value  of  the  nostulate  can  be  determined# 
Stability  in  the  large  is  an  especially  troublesome  point*  Nevertheless, 
the  uniqueness  of  solution  assured  by  the  definition  is  a  good  indication 
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FIGURE  I.  AN  ELASTIC  MATERIAL  (REVERSIBLE):  CURVE 

RISING  TO  THE  RIGHT  Cf6>0 


FIGURE  2.  AN  ELASTIC-PLASTIC  MATERIAL  (TIME 

INDEPENDENT):  RISING  CURVE  0S  >0 
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FIGURE  3.  VISCOELASTIC  MODELS 


(a)  lb) 


FIGURE  4.  EXISTING  SYSTEM  (a)  AND  EXISTING  SYSTEM 

PLUS  EXTERNAL  AGENCY  (b) 
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FIGURE  5-  THE  F UNDAME r  TAL  DEFINITION  REQUIRES 

CONVEXITY  AND  NORMALITY  IN 
PLASTICITY 


FIGURE  6.  PATHS  OF  LOADING-  (Arrows  join  points 

of  the  some  time) 
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FIGURE  7.  [Gj -  C$ 16?  -  &J]  >  0 


FIGURE  8.  (a)  NORMALITY,  6-  =  a^/do^AND  INEQUALITY  [12] 
REQUIRES  CONVEXITY,  lb)  A  CORNER  ALSO  IS  PERMISSIBLE 
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Pasadena  8,  California 
Attn:  Structures  Division  (1) 

Superintendent 

Naval  Post  Graduate  School 

Monterey,  California  (1) 

Commandant 

Marine  Corps  Schools 

Quant ico,  Virginia 

Attn:  Director,  Marine  Corps 

Development  Center  (1) 

Coosrsanding  General 
U.  S.  Air  Force 
Washington  25,  B.  C# 

Attn:  Research  and  Development 

Division  (1) 


Director 

CoBiaander 

Naval  Air  E&cperimental  Station 

Air  Material  Cocssand 

Naval  Air  Material  Center 

Wright -Patterson  Air  Force  Base 

Naval  Base 

Dayton,  Gfrio 

Philadelphia  12,  Pennsylvania 

Attn:  ICREX-3 

Attn:  Materials  Laboratory 

u) 

Structures  Division 

Structures  Laboratory 

u) 

✓ 

Commander 

Officer- in-Charge 

U*  S*  Air  Force  Institute  of 

Underwater  Explosion  Research  Division 

Technology 

Norf  olk  Naval  Shipyard 

Wright-Patterson  Air  Force  Base 

Portsmouth,  Virginia 

Dayton9  Ohio 

Attn:  Dr,  A,  H.  Keil 

Commander 

(2) 

Attn:  Chief,  Applied  Mechanics 
Group 

U.  S.  Naval  Proving  Grounds 

(1) 

Director  of  Intelligence 

Dahlgren,  Virginia 

Headquarters,  U,  S#  Air  Force 
Washington  25*  D*  C# 

Superintendent 

Attn:  P,  V«  Brandi 

Naval  Gun  Factory 

(1) 

(Air  Targets  Division) 

Washington  25,  D.  C« 

(1) 

(1) 


(1) 


(1) 


Commander 

Naval  Ordnance  Test  Station 
Inyokem,  China  Lake,  California 
Attn:  Physics  Division  (1) 

Mechanics  Branch  (1) 

Commanding  Officer  and  Director 
Naval  Engineering  Experiment  Station 
Annapolis,  Maryland  (1) 


Commander 

Air  Force  Office  of  Scientific  Research 
Washington  25,  B#  C. 

Attn:  Mechanics  Division  (1) 

D.  Sm  Atomic  Energy  Connission 
V&shington  25,  B*  G. 

Attn:  Director  of  Research  (2) 
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Director 

National  Bureau  of  Standards 
Washington  25 »  D.  C« 

Attns  Division  of  Mechanics  (1) 

Engineering  Mechanics 
Section  (1) 

Aircraft  Structures  (1) 

Commandant 

U,  S#  Coast  Guard 

1300  E.  Street,  N.W. 

Washington  25 »  D*  C. 

Attns  Chief,  Testing  and 

Development  Division  ( 1 ) 

U#  S#  Maritime  Administration 
General  Administration  Office 
Building 

Washington  25,  D»  C, 

Attns  Chief,  Division  of 

Preliminary  Design  (1) 

National  Advisory  Committee 
for  Aeronautics 
1512  H  Street,  N.W. 

Washington  25,  D.  C. 

Attns  Loads  and  Structures 

Division  (2) 

Director 

Langley  Aeronautical  Laboratory 

Langley  Field,  Virginia 

Attn:  Structures  Division  (2) 

Director 

Forest  Products  Laboratory 
Madison,  Wisconsin  (1) 

Civil  Aeronautics  Administration 
Department  of  Commerce 
Washington  25,  D«  C. 

Attns  Chief,  Aircraft  Engineering 
Division  (1) 

Chief,  Airframe  and  Equipment 
Branch  (l) 

National  Sciences  Foundation 
1520  H  Street,  N«W« 

Washington,  D*  C* 

Attns  Engineering  Sciences  . 

Division  (I) 


National  Academy  of  Sciences 
2101  Constitution  Avenue 
Washington  25,  D*  C» 

Attns  Technical  Director,  Committee 
on  Ships*  Structural  Design  (1) 
Executive  Secretary,  Committee 
on  Undersea  Warfare  (1) 

Professor  Lynn  S*  Beedle 
Fritz  Engineering  Laboratory 
Lehigh  University 

Bethlehem,  Pennsylvania  (1) 

Professor  R*  L.  BispXinghoff 
Department  of  Aeronautical  Engineering 
Massachusetts  Institute  of  Technology 
Cambridge  39,  Massachusetts  (1) 

Professor  Bleich 

Department  of  Civil  Engineering 

Columbia  University 

New  York  27,  New  York  (1) 

Professor  A.  Boiey 

Department  of  Civil  Engineering 

Columbia  University 

New  York  27,  New  York  (1) 

Professor  G#  F*  Carrier 
Pierce  Kali 
Harvard  University 

Cambridge  38 f  Massachusetts  (1) 

Professor  Herbert  Deresiewicz 
Department  of  Civil  Engineering 
Columbia  University 
632  W.  125th  Street 

New  York  27,  N.  Y.  (1) 

Professor  D«  C*  Drucker,  Chairman 
Division  of  Engineering 
Brow^  University 

Providence  12,  Rhode  Island  (l) 

Professor  A*  C«  Eringen 

Department  of  Aeronautical  Engineering 

Purdue  University 

Lafayette ,  Indiana  ( 1 ) 

Professor  W.  Flugge 

Department  of  Mechanical  Engineering 

Stanford  University 

Stanford,  California  (1) 
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Professor  J •  N*  Goodier 
Department  of  Mechanical 
Engineering 
Stanford  University 
Stanford,  California  (1) 

Professor  L*  £*  Goodman 
Engineering  Experiment  Station 
University  of  Minnesota 
Minneapolis  ,  Minnesota  (1) 

Professor  M®  Hetenyi 
The  Technological  Institute 
Northwestern  University 
Evanston,  Illinois  (1) 

Professor  P®  G.  Hodge 
Department  of  aeronautical 
Engineering 

Illinois  Institute  of  Technology 
Chicago  16,  Illinois  (1) 

Professor  M®  J*  Hoff 
Department  of  Aeronautical 
Engineering 
Stanford  University 
Stanford,  California  (1) 

Professor  v.1®  H®  Hoppmann,  II 
Department  of  Mechanics 
Rensselaer  Polytechnic  Institute 
Troy  s  New  York  (1) 

Professor  Bruce  G®  Johnston 

University  of  Michigan 

Ann  Arbor,  Michigan  (1) 

Professor  J®  Kempner 
Department  of  Aeronautical 

Engineering  and  Applied  Mechanics 
Polytechnic  Institute  of  Brooklyn 
99  Livingston  Street 
Brooklyn  2,  New  York  (1) 

Professor  H®  L®  Langhaar 
Department  of  Theoretical  and 
Applied  Mechanics 
University  of  Illinois 
Urbana,  Illinois  (1) 

Professor  B#  J.  Lazan,  Director 
Engineering  Experiment  Station 
University  of  Minnesota 
Minneapolis  14,  Minnesota  (!) 


Professor  £®  H*  Lee,  Chairman 
Division  of  Applied  Mathematics 
Brown  University 

Providence  12,  Rhode  Island  (1) 

Professor  George  H®  Lee 
Director  of  Research 
Rensselaer  Polytechnic  Institute 
Troy,  New  York  (1) 

Mr*  M®  H®  Leiacoe 
Southwest  Research  Institute 
8300  Culebra  Road 

San  Antonio  6,  Texas  (1) 

Professor  Paul  Lieber 

Geology  Department 

University  of  California 

Berkeley  4,  California  (1) 

Professor  Hsu  Lo 
School  of  Engineering 
Purdue  University 

Lafayette,  Indiana  (1) 

Professor  R®  D®  HindTin 

Department  of  Civil  Engineering 

Columbia  University 

632  W®  125th  Street 

New  York  27,  New  York  (1) 

Dr®  A*  Nadai 

136  Cherry  Valley  Road 

Pittsburgh  21,  Pennsylvania  (1) 

Professor  Paul  H®  Naghdi 
Department  of  Engineering 
Mechanics 

University  of  Michigan 

Ann  Arbor,  Michigan  (1) 

Professor  William  Hash 
Department  of  Engineering  Mechanics 
University  of  Florida 
Gainesville,  Florida  (1) 

Professor  E*  H®  Ifewnsark,  Head 
Department  of  Civil  Engineering 
University  of  Illinois 
Urbana,  Illinois  (1) 

Professor  Aris  Phillips 
Department  of  Civil  ^-glnearing 
15  Prospect  Street 
Yale  University 

Now  Haven*  ComBstlmi  (X) 
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Professor  W«  Prager,  Chai,rman 
Physical  Sciences  Council 
Brown  University 

Providence  12,  Rhode  Island  (1) 

Professor  E.  Reissner 
Department  of  Mathematics 
Massachusetts  Institute  of  Technology 
Cambridge  39,  Massachusetts  (1) 

Professor  M,  A*  Sadowsky 
Department  of  Mechanics 
Rensselaer  Polytechnic  Institute 
Troy,  New  York  (l) 

Professor  J.  Stallmeyer 
Department  of  Civil  Engineering 
University  of  Illinois 
Urbana,  Illinois  (1) 

Professor  Eli  Sternberg 
Department  of  Applied  Mathematics 


Brown  University 

Providence  12,  Rhode  Island  (1) 

Professor  S.  P.  Timoshenko 
School  of  Engineering 
Stanford  University 
Stanford,  California  (1) 

Professor  A*  S*  Velestos 

Department  of  Civil  Engineering 

University  of  Illinois 

Urbana,  Illinois  (l) 

Professor  Enrico  Volterra 
Department  of  Engineering  Mechanics 
University  of  Texas 
Austin,  Texas  (1) 

Professor  Dana  Young 
Yale  University 

New  Haven,  Connecticut  (1) 


Dr.  John  F.  Brahtz 
Department  of  Engineering 
University  of  California 
Los  Angeles,  California  (1) 

Mr.  Martin  Goland,  Vice  President 
Southwest  Research  Institute 
3500  Culebra  Road 

San  Antonio,  Texas  (1) 

Mr.  S*  Levy 

General  Electric  Research  Laboratory 

6901  Elmwood  Avenue 

Philadelphia  42,  Pennsylvania  (1) 

Professor  B.  Budiansky 
Department  of  Mechanical  Engineering 
School  of  Applied  Sciences 
Harvard  University 

Cambridge  3$,  Massachusetts  (1) 

Professor  H*  Kolsky 
Division  of  Engineering 
Brown  University 

Providence  12,  Rhode  Island  (1) 

Professor  E*  Orowan 
Department  of  Mechanical  Engineering 
Massachusetts  Institute  of 
Technology 

Cambridge  39,  Massachusetts  (1) 


